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SUMMARY 

Applications  of  computational  aerodynamics  to  aeronautical  research,  design,  and  analysis  have 
increased  rapidly  over  the  past  decade,  and  these  applications  offer  significant  benefits  to  aeroelas- 
ticians.  This  paper  traces  the  past  developments  by  means  of  a  number  of  specific  examples,  and  projects 
the  trends  over  the  next  several  years.  The  crucial  factors  that  limit  the  present  capabilities  for 
unsteady  analyses  are  identified;  they  include  computer  speed  and  memory,  algorithm  and  solution  methods, 
grid  generation,  turbulence  modeling,  vortex  modeling,  data  processing,  and  coupling  of  the  aerodynamic 
and  structural  dynamic  analyses.  The  prospects  for  overcoming  these  limitations  are  presented,  and  many 
improvements  appear  to  be  readily  attainable.  If  so,  a  complete  and  reliable  numerical  simulation  of  the 
unsteady,  transonic  viscous  flow  around  a  realistic  fighter  aircraft  configuration  could  become  possible 
within  the  next  decade.  The  possibilities  p*  using  artificial  intelliience  concepts  to  hasten  the  achieve¬ 
ment  of  this  goal  are  also  discussed. 

1.  INTRODUCTION 

The  extraordinary  growth  in  computer  technology  of  the  past  two  decades  has  revol utioni ted  the  design 
if  modern  aircraft,  especially  aircraft  that  cruise  or  maneuver  in  the  transonic  flow  '■egime.  There  is  no 
need  to  belabor  the  contributions  of  computational  fluid  dynamics  .CFO)  to  solving  aeronautical  problems, 
as  a  large  number  of  review  papers  nave  been  written  on  this  subject;  Refs.  1-1  are  but  a  small  sample. 
Suffice  It  to  say  that  the  requirements  and  potential  benefits  for  predicting  the  performance  and  steady 
airloads  of  advanced  aircraft  have  been  a  major  driving  force  in  the  development  of  existing  CFO  technology 
and  contemporary  supercomputers;  and  that,  as  a  result,  impressive  capabilities  exist  today.  This  paper 
addresses  the  current  status  and  the  prospects  for  developing  new  CFO  methodology  to  predict  the  aeroelas- 
tic  behavior  of  future  advanced  aircraft.  <Je  will  focus  on  the  principal  factors  that  will  determine  the 
success  or  failure  of  computational  aerodynamics  to  meet  the  needs  of  aeroelasticians ,  and  on  the  future 
developments  that  might  be  expected  to  alter  these  factors. 

The  demands  of  the  aeroelasticity  community  on  aerodynamici sts  are  staggering.  For  dozens  or  even 
nundreds  of  cases,  they  are  asked  to  provide  the  three-dimensional,  unsteady  airloads  on  complex  geometries 
including  external  stores)  at  transonic  flight  conditions  for  which  viscous  effects  are  important,  and  to 
couple  these  aerodynamic  characteristics  with  the  complex  structural-dynamic  behavior  of  the  airframe. 
Furthermore,  these  results  are  to  be  obtained  ouickly,  and  without  significantly  increasing  the  national 
debt.  Thus,  calculating  the  flutter  boundaries  of  flexible  aircraft  flying  at  transonic  speeds  is  at  least 
an  order  of  magnitude  more  difficult  than  the  more  widely  publicized  task  of  predicting  the  performance  of 
aerodynamic  configurations  in  steady  flight.  Finally,  aeroelastic  calculations  are  more  specialized,  and 
"ewer  people  are  working  in  this  area,  than  in  steady  aerodynamics.  Fortunately,  there  is  clear  evidence 
•n  the  16  papers  of  this  Specialists'  Meeting,  and  elsewhere,  that  considerable  progress  has  been  and  is 
being  made  toward  adapting  and  extending  the  latest  CFO  methodology  for  steady  flows  into  the  unsteady 
domain.  However,  we  will  demonstrate  that  new  thrusts,  new  ideas,  and  new  levels  of  effort  will  be 
'•eguired  to  meet  the  aeroelasticians'  future  requirements. 

One  way  to  view  the  task  of  determining  flutter  behavior  is  to  think  of  it  as  an  exercise  in  simula¬ 
tion.  That  is,  we  want  to  avoid  unacceptable  catastrophes  in  real  life  by  analyzing  appropriate  risk-free 
simulations  of  the  phenomena.  One  may  note,  in  passing,  that  this  concept  is  more  akin  to  the  issues 
facing  the  nuclear  power  industry  than  to  those  of  predicting  aircraft  performance.  The  point  is  that  we 
wisn  to  simulate  a  complex  dynamics  problem  in  solid  mechanics  which  is  driven  by  complicated,  boundary- 
dependent,  and  nonlinear  aerodynamic  forcing  functions.  We  now  have  available  to  us  highly  developed 
techniques  for  physical  simulations  (which  we  call  wind-tunnel  testing)  and  rapidly  developing  techniques 
For  numerical  simulation  (which  can  combine  large  scientific  computers,  structural -dynamic  analyses,  and 
computational  aerodynamics).  Each  type  of  simulation  has  its  strengths  and  its  limitations,  as  indicated 
'r  fiq.  1.  The  relative  merits  are  frequently  .debated;  but  we  shall  not  dwell  on  these  issues  here,  except 
to  argue  for  the  judicious  use  both  of  experiment  and  computation  to  complement  each  other. 

As  discussed  in  Section  3,  our  projection  is  that  complete  and  reliable  numerical  simulations  will 
become  possible  within  the  next  decade  for  complex  configurations,  and  that  high-quality  physical  experi¬ 
ments  will  play  crucial  roles  in  developing  and  validating  these  numerical  simulations.  However,  the  cost 
if  the  complete  simulations  may  well  be  excessive  for  the  hundreds  of  combinations  of  flow  parameters, 
structural  Frequencies  and  mode  shapes,  and  wing-store  configurations  that  aeroelasticians  will  probably 
want  to  analyze.  Therefore,  there  will  clear'y  be  an  ongoing  need  for  less  costly,  more  approximate  simu- 
’ation  methods,  even  though  some  accuracy  may  be  sacrificed.  The  development  and  validation  of  such  engi¬ 
neering  simulations  will  be  greatly  enhanced  by  an  intelligent  combination  of  the  large-scale  numerical 
and  physical  simulations. 

The  main  thrusts  of  this  paper,  then,  are  to  demonstrate  the  rapid  and  continuing  growth  of  computa¬ 
tional  aerodynamics,  to  indicate  the  principal  areas  that  must  be  further  developed  if  computational  aero- 
dynamicists  are  to  provide  significantly  better  tools  for  aeroelasticians,  and  to  examine  the  requirements 
for  a  complete,  time-accurate  numerical  simulation  of  the  unsteady,  transonic  viscous  flows  around  a 
■■ealistic  Fighter  aircraft  configuration.  This  paper  should  be  considered  to  be  complementary  to  the 
broader,  companion  review  of  V.  L.  Peterson,  entitled  'Trends  ’n  Computational  Capabilities  for  Fluid 
Dynamics  '  (Ref.  3),  given  earlier  at  this  Specialists'  Meeting.  No  attempt  has  been  made  to  review  com¬ 
prehensively  the  large  and  rapidly  growing  body  of  literature  on  unsteady  computational  aerodynamics. 
Rather,  our  somewhat  random  and  parochial  choice  of  representative  examples  largely  reflects  the  research 
with  which  we  are  most  familiar  and  the  results  which  are  most  readily  available  to  us. 
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Fig.  1.  Unsteady  flow  reginws  for  transonic  wings. 

at  lower  reduced  frequencies,  where  the  shock-wave  motion  is  the  largest  and  the  resultant  unsteady  aspects 
of  the  fluid  physics  are  the  most  complex. 

Unfortunately,  the  current  computational -aerodynamics  codes  that  might  be  capable  of  capturing  these 
low-frequency  complexities  tend  to  have  severe  stability  limitations  with  respect  to  the  maximum  time  step 
that  can  be  used.  This  translates  into  long  computational  times  per  cycle  of  oscillation.  In  addition  to 
purely  financial  reasons,  long  CPU  time  impedes  progress  in  several  ways.  First,  it  discourages  new  users 
from  trying  unfamiliar  codes  and  becoming  comfortable  with  them;  second,  it  discourages  experienced  users 
from  experimenting  with  the  codes  in  new  applications;  and  third,  it  limits  code  developers  in  their  efforts 
to  refine,  tune,  and  extend  the  methods. 

An  impractical  conclusion  that  could  be  drawn  from  Fig.  2  is  that  aeroelastic  ai  ilyses  would  be  simpli¬ 
fied  if  the  vehicles  were  restricted  to  subsonic  speeds,  or  if  the  natural  frequencies  of  the  structures 
could  be  increased  by  an  order  of  magnitude. 

2.1  Representative  Calculations-  1974  and  1984 

A  better  understanding  of  the  current  trends  in  computational  aerodynamics  can  be  obtained  by  tracing 
the  growth  in  capabilities  that  has  occurred  over  the  past  decade.  In  this  section  we  will  note  a  few 
examples  that  were  particularly  noteworthy  as  state-of-the-art  circa  1974,  and  discuss  the  corresponding 
capabilities  today. 

2.1.1  Steady  flow,  complex  geometries 

For  many  years,  linear  panel  methods  have  been  the  primary  tool  for  analyzing  complete  vehicles  of 
complex  geometry.  Figure  3  illustrates  the  degree  of  sophistication  that  has  been  and  is  possible.  Larger 
computers  have  enabled  the  use  of  more  panels,  with  a  corresponding  improvement  in  the  resolution  of  the 
surface  airloads.  However,  the  primary  advance  in  panel  methods  since  1974  is  the  present  capability  to 
treat  supersonic  problems  with  the  same  surface  representation  that  before  could  be  done  only  for  subsonic 
cases.  Figure  4  shows  representative  results  (Ref.  9)  for  a  fighter  configuration  with  canards.  The  lift 
of  the  aircraft  is  well  predicted,  but  further  improvements  are  required  for  predicting  the  drag. 

Insofar  as  nonlinear  methods  are  concerned,  three-dimensional  transonic  smal 1 -disturbance  calculations 
of  wing-body  comb'nations  were  possible  with  the  Bai ley-Bal lhaus  code  (Ref.  10)  a  decade  ago.  More  compli¬ 
cated  configurations  are  routinely  analyzed  with  full  potential  methods  today,  and  Euler  methods  are  coming 
to  the  forefront;  e.g..  Refs.  11-13.  Figure  5  shows  the  recent  calculations  of  Jameson  and  Baker  (Ref.  11). 
These  results  have  not  been  verified  by  independent  calculations  or  comparison  with  experiment,  and  the 
very  coarse  grid  on  the  tail  surfaces  is  probably  inadequate  for  resolving  the  flow  in  that  region.  How¬ 
ever,  the  influence  of  the  body  and  tail  on  the  flow  over  the  wing  is  probably  captured  accurately  enough 
'n  this  simulation. 

2.1.2  Complex  steady  flow,  simple  aerodynamic  shapes 

Jameson's  FL06  transonic  potential-flow  code  (Ref.  14)  for  airfoils  with  shock  waves  came  into  general 
use  in  the  early  1970s;  and  by  1974,  weak  viscous  corrections  had  been  added  (Ref.  15).  With  regard  to 
viscous-dominated  flows,  in  that  era  Mehta  (Ref.  16)  treated  the  fully  separated  flow  of  an  airfoil  at  high 
angle  of  attack  using  the  laminar  Navier-Stokes  equations,  and  computed  the  self-induced  fluctuations  as 
well  as  the  mean  airloads.  Today  the  stalled  airfoil  at  high  Reynolds  numbers  remains  an  unsolved  problem, 
but  this  IS  mainly  because  ■^f  the  turbulence  modeling,  and  not  because  of  the  computational  barriers. 

A  decade  ago  Oeiwert  (Ref.  17)  treated  shock- induced  separation  on  a  nonliftinq  airfoil  with  the 
Reynolds-averaged  Navier-Stokes  equations,  using  an  algebraic  mi «ing-lenqth  model  of  the  turbulence.  That 
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capability  nas  since  been  extended  to  a  lifting  transonic  wing  by  Mansour  (Ref.  18),  including  both  shock 
waves  and  tip  vortex  formation  (Fig.  6),  and  to  afterbodies  with  propulsive  jets  by  Deiwert  and  Rothmund 
(Ref.  19)  (Fig.  7). 
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-■g.  ’  ’hin-layer,  Reynolds-averaged  Navier-Stokes  calculations  of  a  conical  afterbody  with  a  propulsive 
;et  Ref.  19);  •  2.0,  x  »  6",  =  2.5.  Re  =  1.5  «  10®,  140  »  100  x  20  grid. 

Although  Mansour's  wing  calculations  were  performed  on  a  relatively  coarse  grid,  several  hours  of  CPU 
time  were  required  to  obtain  the  results  shown  in  Fig.  6.  Furthermore,  the  agreement  with  the  experimental 
data  IS  only  fair,  and  probably  10  times  as  many  grid  points  would  be  required  to  resolve  the  details  of 
the  flow.  Nevertheless,  this  investigation  represents  a  milestone  in  analyzing  viscous  wing  flows. 

The  afterbody  calculations  shown  in  Fig.  7  provide  a  remarkable  amount  of  detail  by  using  only  a 
modest  number  of  grid  points.  This  capability  stems  from  a  high  degree  of  specialized  experience  with 
this  particular  code  and  this  particular  class  of  problems,  and  it  illustrates  the  value  of  having  skilled 
experts  to  work  with  a  family  of  codes.  Even  so.  difficult  cases  with  larger  separation  zones  still  give 
problems  (Ref.  20). 

2.1.3  Unsteady  inviscid  flow 

One  of  the  showcase  results  of  the  mid-1970s  was  the  calculation  of  Magnus  and  Yoshihara  (Ref.  21) 
for  an  oscillating  airfoil  with  a  strong  shock  wave,  using  an  explicit  Euler  method  (Fig.  8).  Another 
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Fig.  3.  Euler  calculations  of  a  transonic  oscillating  airfoil  (Ref.  21); 

NACA  64A-410  airfoil.  M  *  0.72,  X/U,  =  0.20. 

oioneering  set  of  results  was  obtained  by  Caradonna  and  Isom  (Ref.  22),  using  a  three-dimensional,  unsteady 
transonic  smal 1 -disturbance  code,  shown  in  Fig.  9.  The  two-dimensional  small -disturbance  code  LTRAN2 
(Ref.  23)  also  became  available  soon  afterward,  and  it  has  been  used  extensively  ever  since.  Each  of 
these  helped  pave  the  way  for  the  computational  capabilities  that  exist  today,  and  they  provided  insights 
into  unsteady  effects  that  could  not  have  been  obtained  either  by  linear  theory  or  by  experiments. 
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Fig.  9.  Transonic  smal 1 -disturbance  calculations  of  a  nonlifting  helicopter  blade  tip  (Ref.  22). 


The  unsteady  transonic  smal 1 -disturbance  method  has  since  achieved  a  high  level  of  maturity  in  both 
two  and  three  dimensions.  Also,  approximate  viscous  corrections  have  been  added,  and  the  aerodynamic  cal¬ 
culations  have  been  coupled  with  the  structure  (Refs.  24-27).  Figure  10  shows  the  results  of  Guruswamy  and 
Goorjian  (Ref.  27)  for  a  low-aspect-ratio  oscillating  wing.  These  calculations,  using  51,200  grid  points 
and  1024  time  steps  per  cycle,  required  about  30  min  of  CPU  time  on  a  Cray  XMP  computer  to  compute  three 
cycles  of  oscillation. 
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Fig.  lO.  Transonic  small-disturbance  calculations  of  the  oscillating  F-5  wing  (Ref.  27); 

M  =  0.90.  K  =  -c/U.  «  0.55. 

The  ful 1 -potential  results  of  Malone  et  al.  (Ref.  28)  for  the  same  oscillating  wing  are  shown  in 
Fig.  11.  In  this  case,  only  18,000  grid  points  were  used,  with  correspondingly  less  resolution.  The  cal¬ 
culations  reduired  several  hours  on  a  VAX  11/780  computer.  Finally,  the  ful 1 -potential  calculations  by 
Isogai  and  Suetsugu  (Ref.  29).  for  wings  with  part-span  oscillating  flaps,  may  also  be  mentioned  to  illus¬ 
trate  the  current  capabilities  for  inviscid  flows. 

2.1.4  Unsteady  viscous  flow 

Although  unsteady  effects  on  turbulent  boundary  layers  were  computed  by  several  investigators  in  the 
early  1970s,  there  were  essentially  no  unsteady  viscous  results  available  10  years  ago  that  were  of  direct 
interest  to  aeroelasticians.  Today  we  can  point  to  studies  of  oscillating  airfoils  and  flaps  using  the 
Reynolds-averaged  Navier-Stokes  equations  with  simple  eddy-viscosity  turbulence  models  that  can  be  run  in 
an  hour  or  less  on  modern  supercomputers.  As  an  example,  the  calculations  of  Horiuti  et  al.  (Ref.  30)  are 
shown  in  Fig.  12.  This  investigation  also  includes  a  study  of  the  effects  of  wind-tunnel  walls,  which  can 
significantly  alter  the  phase  of  the  unsteady  pressure  distribution  behind  the  shock  wave. 

Another  recent  Navier-Stokes  calculation  at  transonic  speed  and  a  high  Reynolds  number  (Fig.  13), 
from  Ishii  and  Kuwahara  (Ref.  31),  illustrates  the  growing  CFD  capability  in  Japan.  It  is  interesting  to 
note  that  these  results,  which  have  not  been  validated  by  experimental  comparisons,  were  obtained  with  no 
turbulence  model  at  all.  Finally,  the  recent  review  of  compressible  Navier-Stokes  solutions  by  Shang 
(Ref.  32)  may  be  consulted  for  further  examples  and  for  a  comprehensive  bibliography. 

2.2  Cost  and  Capability  Trends 

The  preceding  examples  indicate  the  growing  capabilities  to  solve  challenging  aerodynamics  problems. 
Much  of  this  progress  can  be  traced  directly  to  the  extraordinary  growth  in  computer  technology,  as  dis¬ 
cussed  in  Refs.  1-8.  Computer  speed,  memory  size,  and  cost  are  all  important  factors  in  assessing  the 
present  and  future  capabilities  for  performing  complex  aeroelastic  analyses.  An  overview  of  the  trends 
for  these  factors  is  given  in  Figs.  14-16,  from  Peterson  (Ref.  8). 

Figures  14  and  15  show  that  speed  and  memory  capacity  continue  to  grow  more  rapidly  than  the  costs  of 
the  machines.  Consequently,  the  relative  cost  of  performing  aerodynamic  calculations  is  decreasing  dramati¬ 
cally  (Fig.  16).  The  improvements  in  algorithms  and  methods  of  analysis  are  more  difficult  to  quantify  than 
are  those  in  hardware;  but  the  general  trends  are  clear,  and  researchers  are  confident  of  further  gains. 

The  net  result  is  that  the  cost  of  performing  a  given  computation  has  decreased  three  or  more  orders  of 
magnitude  per  decade  (Refs.  6-8),  and  this  trend  is  projected  to  continue  for  some  time. 
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Fig.  n.  Unsteady  full-potential  calculations  of  the  oscillating  F-5  wing  (Ref.  28); 

M  =■  0.90,  K  *  „c/U.  •  0.55. 

An  even  irore  important  trend,  which  is  difficult  to  portray  graphically,  is  the  strong  tendency  to 
attack  increasingly  difficult  problems  {and  with  greater  confidence)  as  speed,  memory,  numerical  methods, 
and  physical  modeling  improve.  To  paraphrase  the  oral  version  of  Ref.  3,  concerning  the  impact  of  CFO  on 
conFtercial  aircraft  design. 

The  total  costs  of  computing  have  gone  up,  not  down-we  just  do  a  lot  more  CFD  than  ever  before 

in  the  quest  for  superior,  innovative  designs. 

This  trend  is  also  implicit  in  MacCormack’s  prediction  (Ref.  33)  that  a  Reynolds-averaged  Navier-Stokes 
solution  for  a  complete  aircraft  will  be  obtained  in  1985. 

2.3  Summary  of  Current  Status 

The  growth  of  CFD  over  the  past  decade  has  given  us  adequate  capabilities  to  model  compressible  flows 
with  imbedded  shock  waves  and  weak  viscous  effects  (i.e.,  at  low  angle  of  attack  and  without  shock-induced 
separation).  This  can  be  done  on  simple,  but  practical,  bodies  undergoing  smal 1 -ampl i tude  motions,  in 
free  air  or  solid-wall  wind  tunnels.  However,  many  of  the  codes  have  not  been  adequately  validated  and 
calibrated.  Also,  we  still  have,  at  best,  only  marginal  capabilities  for  strong  vortices,  strong  turbulent 
viscous  effects,  complex  geometries  in  the  transonic  regime,  or  simulation  of  venti lated-wal 1  wind  tunnels. 

The  combined  hardware  and  software  costs  of  computing  today's  problems  are  not  trivial,  especially 
when  user-manpower  costs  are  included.  However,  the  relative  costs  of  computations  have  dropped  steadily 
by  a  factor  of  about  1000  per  decade  over  the  past  20  years.  In  addition,  the  growth  in  capability  and 
the  reduction  in  turn-around  time  for  a  typical  calculation  are  even  more  important  than  the  cost  trends 
to  many  segments  of  the  aircraft  industry.  Except  for  the  limitations  of  turbulence  modeling,  the  show¬ 
case  problems  of  1974  can  be  solved  routinely  in  1984. 


VELOCITY  VECTOBS  DENSITY  CONTOURS 


Compressible  flow  around  a  circular  cylinder  at  M  =  0.8  and  8e  =  5  •  10'’  ;Ref.  31 
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Fig.  14.  The  growth  of  computer  soeed  and  cost  (Ref.  8). 
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Fig.  15.  The  growth  of  computer  memory  (Ref.  8). 
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3.  THE  CHALLENGES  FOR  AEROELASTIC  APPLICATIONS 


Despite  the  progress  that  has  been  made  in  computational  aerodynamics,  there  are  definite  limits  to 
what  can  be  done  practically,  especially  in  aeroelastic  applications.  Also,  the  computing  cost  for  the 
aerodynamics  part  of  a  typical  flutter  simulation  for  a  complex  vehicle  flying  at  transonic  speeds  could 
exceed  the  structural-dynamics  part  by  a  considerable  margin.  In  this  section,  we  turn  to  some  of  the 
pacing  items  and  current  limitations  of  CFO  as  they  relate  to  aeroelastic  problems.  He  shall  also  consider 
the  extent  to  which  some  of  these  limitations  might  be  relaxed  or  eliminated,  and  indicate  some  areas  in 
which  new  developments  could  pay  rich  dividends. 

3.1  Limiting  Factors  for  Computational  Aerodynamics 

Figure  17  illustrates  the  most  important  pacing  items  in  applying  computational  aerodynamics  to  aero¬ 
elastic  problems.  Grid  generation,  turbulence  modeling,  and  computer  hardware  and  software  were  identifieo 
in  Refs.  5  and  6  as  primary  pacing  items,  and  this  will  remain  true  for  the  foreseeable  future  in  almost 
all  areas  of  CFD.  Current  applications  are  limited  by  both  computer  speed  and  memory;  the  relative  impor¬ 
tance  of  each  can  be  debated  for  aircraft  performance  predictions  and  design.  However,  as  we  shall  see, 
the  roles  of  the  algorithm  and  solution  methodology  tend  to  be  more  important  in  time-dependent  problems 
than  in  steady-flow  predictions.  This  is  because  smaller  time  steps  are  often  needed  for  stability  con¬ 
siderations  than  for  accurately  capturing  the  unsteady  features  of  the  flow. 


GRID  GENERATION  TURBULENCE  MODELING 


COMPUTER  POWER 


AEROELASTICITY  -  LACK  OF  CFD-SKILLED  ENGINEERS  AND  MANAGERS 
Fig.  17.  Pacing  items  in  computational  fluid  dynamics. 


Another  novel  aspect  of  aeroelasticity  is  less  of  a  technical  issue  than  a  management  one;  namely, 
the  shortage  of  engineers  and  research  scientists  who  are  specialists  in  both  computational  'luid  dynamics 
and  structural  dynamics.  To  a  certain  extent,  the  rate  of  progress  can  be  expected  to  be  proportional  to 
the  level  of  effort  expended  and  to  the  skills  of  the  people  exerting  the  effort.  ;n  addition,  there  are 
even  fewer  managers  who  have  been  trained  in  both  these  disciplines. 

Insights  into  the  magnitude  of  the  challenges,  and  prospects  ^or  nonlinear  aeroelastic  applications, 
can  oe  obtained  by  analyzing  the  factors  that  determine  the  solution  times,  or  CPU  nequi  rements ,  ''or  a  typi¬ 
cal  time-accurate,  unsteady  airloads  calculation.  For  most  CFD  methods,  the  CPU  time  can  be  expressed  as 

CPU  =  A  <  -  1^,  FLOPS 

where 

A  =  "numerical  inefficiency’  factor 

=  number  of  operations  per  grid  point  oer  time  step 
Np  =  number  of  grid  points 

m  =  1  for  finite-difference,  *or  panel  methods 

Nj  =  number  of  time  steps  =  (number  of  time  steos.'c/c le '  <  I'number  of  cycles!  =  number  ,-eference 
lengths/cycle)  '  (number  of  cycles'/'.; 

!  =  nondimensional  time  step  =  "J  .'.t/L 
FLOPS  »  number  of  floating-point  arithmetic  operations  per  unit  time 


Here  we  introduce  the  efficiency  factor  A,  to  emphasize  that  the  code  may  not  take  full  advantage  of 
the  computer  being  used;  in  practice  it  is  a  function  of  the  programning  efficiency,  the  degree  of  vectori- 
zation,  the  coupling  between  the  grid  and  the  solution  algorithm,  the  user  experience,  etc.  Ideally,  its 
value  should  approach  unity;  but  especially  with  the  advent  of  supercomputers  with  novel  architecture,  it 
could  be  much  larger. 

The  number  of  arithmetic  operations  per  grid  point  per  time  step,  Wqj,  is  a  strong  function  of  the 
numerical  method;  that  is,  of  the  flow  equations,  the  boundary  conditions,  the  solution  algorithm,  and  the 
grid.  The  quantity  Nr  represents  the  number  of  grid  points  for  a  finite-difference  method,  the  number 
jf  elements  for  a  finite-element  method,  or  the  number  of  panels  for  a  panel  method.  Consequently,  WqjNq 
represents  the  number  of  arithmetic  operations  that  must  be  performed  at  each  Iteration  or  time  step; 
although,  in  some  instances  with  panel  methods,  Nq  log  Nq  is  a  more  accurate  representation  than  N^. 

Ideally  for  aeroelastic  applications,  the  total  number  of  time  steps,  Nj,  would  simply  be  the  number 
of  time  steps  per  cycle  multiplied  by  the  number  of  cycles  needed  to  determine  the  flutter  characteristics, 
-owever,  many  nonlinear  aerodynamics  codes  have  stability  or  accuracy  limits  that  are  determined  by  a  non- 
limensional  time  step,  ct  =  U  ct/L.  Thus  the  maximun  permissible  value  of  or  typically  depends  upon  the 
.o’-r.ievity  of  the  problem,  the  algorithm,  the  grid,  and  the  desired  accuracy. 

-■nally,  the  computing  speed,  FLOPS,  is  a  function  of  the  computer  clock  speed  and  architecture,  the 
iata  nanagement  techniques  of  the  code,  the  memory  requirements  (in-core  or  external  memory),  and  the  solu¬ 
tion  algorithm,  ’'hus  it  is  clear  that  many  different  factors  determine  the  CPU  time,  and  the  cost,  of  an 
ierodynamic  calculation. 

estimates  nave  been  made  of  the  solution  times  that  would  be  required  to  run  a  wide  range  of  contem- 
corary  time-accurate  methods  on  modern  supercomputers.  Table  2  shows  a  breakdown  of  the  factors  in  Eq.  (1) 
•-.r  a  winq  of  moderate  complexity  undergoing  three  cycles  of  oscillatory  motion  and  25  chord  lengths  of 
•.rave’  per  cycle,  running  on  a  computer  with  a  nominal  sustained  rate  of  80  million  floating-point  opera- 
•'.-''s  per  second.  By  moderate  complexity,'  we  mean  something  a  bit  more  complicated  than  the  wing  shown 
r  0,  This  might  include,  for  example,  a  relatively  clean  wing-body  combination,  a  wing  with  a  flap, 
a  wing  with  a  tip  tank  or  tip-mounted  rocket,  etc.,  but  not  a  wing  with  multiple  external  stores.  Other 
assumptions  are  noted  in  the  notes  to  Table  2. 


'ABLE  2.  COMPUTATIONAL  REQUIREMENTS  FOR  COMPLEX  OSCILLATING  WINGS 


Flow  model 

'“^GT 

'‘g 

pt 

CPU, 

minutes 

Memory, 
million  words 

Notes 

Nonl inear  panel 

‘^GT  ■ 

s 

2  X  10“ 

0.05 

60 

2.0 

* 

Small  disturbance 

100 

10‘ 

0.06 

8 

0.6 

Full  potential 

600 

10* 

0.04 

23 

2.0 

Full  potential 

and  integral  B.L. 

630 

10* 

0.02 

50 

2.0 

Euler 

3000 

10* 

0.01 

450 

3.0 

i,v 

Euler  and  finite 

difference  B.L. 

2000 

2 

.  10* 

0.01 

500 

5.0 

Thin-layer 

Navier-Stokes 

3600 

10‘ 

0.005 

1 1  .000 

30 

Reynolds-averaged 

Navier-Stokes 

4500 

2 

.  10* 

0.004 

35,000 

60 

Notes:  4*2 

i.  A  =  3 

at  for  time  accuracy 
i.  pt  for  stability  limitations 

Wqt  includes  100  for  grid  generation 
Sw  increase  in  Wqj  for  boundary  layer 
y.  Wqj  =  SCO  in  viscous  layer,  3000  in  inviscid  region 

It  should  be  mentioned  that  a  time  step  limit  of  at  =  0.05  has  been  assumed  as  a  rather  subjective 
estimate  of  what  is  required  to  resolve  accurately  unsteady  transonic  effects,  including  significant  shock 
wave  motion.  For  cases  which  can  be  considered  as  almost  linear  perturbations  about  a  nonlinear  mean  flow, 
much  larger  values  might  suffice.  We  also  acknowledge  that  our  estimates  of  the  stability  limitations  on 
pt  are  very  approximate,  and  the  numbers  given  are  intended  to  give  a  sense  more  of  the  relative  values  of 
the  various  methods  than  of  the  absolute  values.  The  important  point  is  that,  currently,  the  more  sophis¬ 
ticated  the  method,  the  more  severe  is  the  stability  i-estriction  on  pt  for  highly  nonlinear  problems. 

The  sol ut ion- time  requirements  are  compared  graphically  in  Fig.  18.  'We  note  again  that  all  of  these 
results  are  very  approximate,  accurate  to  one  significant  figure  at  best.  In  Fig.  18,  the  "time-1 ineari zed" 
estimate  refers  to  any  of  the  numerous  methods  that  are  available  for  obtaining  the  harmonic  components  of 
the  unsteady  airloads  as  a  linear  perturbation  about  a  nonlinear  mean-flow  condition;  i.e..  Regime  3  in 
Fig.  2.  The  estimates  for  the  nonlinear  panel  methods  (which  are  based  on  Refs.  34-36,  and  private  con¬ 
versations  with  Forrester  .Johnson  of  the  Boeing  Military  Airplane  Company  and  Larry  Erickson  of  NASA  Ames) 
are  even  more  approximate,  as  these  methods  have  not  yet  been  used  to  calculate  the  time-accurate  evolution 
of  transonic  flow  fields.  ,  , 


Fig.  18.  Solution  times  for  a  wing  of  moderate  complexity;  75  chords  travel/case,  1984  algorithms, 

80  mflops. 


Similarly,  the  estimates  for  combining  an  Euler  method  and  finite-difference  boundary  layer  are  purely 
speculative.  They  are  based  on  the  work  of  Van  Dalsem  and  Steger  (Ref.  37)  with  the  steady  full -potential 
equation,  and  private  conversations  with  them;  however,  this  approach  has  yet  to  be  implemented  In  an 
actual  unsteady,  three-dimensional  code.  The  basic  concept  is  very  attractive,  however.  The  additional 
number  of  grid  points  In  the  viscous  layer  and  wake  is  assumed  to  be  the  same  as  the  number  In  the  Invlscid 

field,  but  the  number  of  arithmetic  operations  per  viscous  grid  point  is  only  about  600.  Therefore,  the 

average  value  of  Wgf  decreases  for  the  coupled  system  of  equations,  and  the  product  of  Wqt  eod  Ng 
increases  only  about  30%. 

The  information  in  Table  2  and  Fig.  18  Implies  that  if  the  transonic  effects  are  mild  enough  to  fall 
within  the  scope  of  the  time- linearized  methods,  then  extensive  flutter  analyses  are  probably  affordable 
now.  But  time-accuracy  is  expensive  today  for  any  nonlinear  method,  and  is  probably  prohibitively  so  for 

anything  more  sophisticated  than  the  potential  methods.  This  raises  the  question  of  what  can  be  done,  or 

what  Is  likely  to  be  done,  over  the  next  few  years  to  bring  the  costs  of  all  the  methods  down  to  acceptable 
levels. 

3.2  Targets  of  Opportunity 

Independent  of  the  efforts  of  the  CFO  and  aeroelasticl ty  communities,  supercomputer  technology  can  be 
expected  to  increase  the  computation  speeds  from  80  to  250  mflops  or  more,  within  the  next  3  years,  and  to 

provide  adequate  memory  to  meet  the  requirements  listed  in  Table  2.  In  addition,  a  number  of  specific 

improvements  can  be  readily  foreseen  in  computational  aerodynamics.  Although  many  of  these  improvements 
will  take  significant  amounts  of  time  and  effort  to  develop  and  to  validate,  they  represent  advances  that 
are  within  our  grasp. 

3.2.1  Algorithm  and  numerical  technique  Improvements 

Equation  (1)  provides  the  basis  for  assessing  the  possibilities  of  significantly  reducing  the  CPU  times 
for  the  various  time-accurate  methods.  The  main  thrusts  should  obviously  be  to  reduce  the  total  number  of 

floating-point  operations  for  each  time  step,  as  indicated  in  Table  3,  and  to  decrease  the  total  number  of 

time  steps  by  increasing  pt  (Table  4).  In  the  following  discussion,  it  will  be  assumed  that  the  previous 
upper  limit  of  Pt  =  0.05,  based  on  accuracy  considerations,  can  be  increased  to  U.IO  by  a  combination  of 
experience  and  improved  grid  techniques. 

Nonlinear  panel  methods  -  These  methods  have  been  under  development  for  steady  flows,  where  multi-grid 
techniques  are  useful  (Refs.  34  and  36).  They  would  seem  to  be  candidates  for  significant  reductions  in 
the  number  of  iterations  required  per  timg  step,  and  for  reducing  the  dependence  on  the  number  of  surface 
panels  or  field  grid  points  to  perhaps  or  Nq  log  instead  of  N‘.  Also,  it  should  be  possible  to 

increase  the  time  step  limit  to  0.1,  as  noted  above,  and  to  remove  the  stability  restrictions  that  seem  to 
affect  current  integral  boundary-layer  treatments  of  viscous  effects. 

Transonic  small  disturbance  methods -On  the  other  hand,  this  approach  has  already  matured  to  the  point 
where  the  number  of  arithmetic  operations  required  is  not  likely  to  decrease  significantly  below  current 
levels.  Room  for  improvement  exists  primarily  in  decreasing  A  to  unity  by  rewriting  the  existing  codes 
'Refs.  24  and  26),  and  by  increasing  .st,  as  above.  This  can  probably  be  done  by  treating  more  of  the 
'secondary'  terms  implicitly  instead  of  explicitly.  With  additional  experience,  the  same  accuracy  can 
probably  be  attained  with  a  30-50%  x-eduction  in  the  total  number  of  grid  points,  with  a  corresponding 
reduction  in  CPU  time. 


TABLE  J.  REDUCING  COMPUTATION  TIMES 
A  ■  “  effective  total  operations/time  step 


Flow  model 

Today 

ops/time  step 

Future 

Nonl inear  panel 

2  «  10* 

Reduce 

m 

“gt 

to  1  iteration/time  step 

Small  disturbance 

3  '  10'  Reduce 

A  to  1 .0 

Full  potential 

6  <10'  Reduce 

'‘(5 

by  50%  with  grid  adaption 

Euler  implicit 

3  <  10*  Reduce 

•<GT 

by  algorithm 

Euler  implicit  plus  boundary 

1  ayer 

A  <  10*  Reduce 

by  grid  adaption 

Euler  expl icit 

0.3  «  10*  Reduce 

by  grid  adaption 

Thin-layer  Navier-Stokes 

4  ’<  10*  Reduce 

■i^GT 

by  algorithm 

Reynolds-averaged  Navier-Stokes 

9  »  10*  Reduce 

by  grid  adaption  and  zonal  modeling 

TABLE 

4.  REDUCING  COMPUTATION 

TIMES 

N,/A- 

=  total  number  of  time  steps 

Flow  solver 

At  today 

Future  improvements 

Nonl i near  panel 

Small  Disturbance 

Full  potential 

0.05 

0.02 

inviscid  Increase  to 

B.L. 

0.1  by  experience  and  grid  adaption 

Eqler  imc’icit  (includes  8.L.' 

0.02 

(stability)  O.IO 

remove  stability  limit) 

Eul  l't  e.vp  1  ici  t 

0.01 

(stability)  (increase  ax)  _ ^ 

'iavier-Stakes  implicit 

'0.0105 

(stability)  0.1  (remove 

stability  limit) 

Full  potential  methods -The  two  ma’n  ways  to  improve  this  approach  are  to  reduce  the  number  of  grid 
points  required  by  using  solution-adaptive  grids,  without  increasing  NgT.  *''<1  to  increase  at  by  means  of 
the  better  grids.  A  likely  additional  improvement  will  be  to  include  an  integral  unsteady  boundary  layer 
formulation,  with  negligible  increase  in  Hqt  And  with  no  stability  restriction  on  at. 

Euler  methods  -  It  may  be  possible  to  reduce  the  value  of  Wqt  by  about  a  third  for  this  approach,  and 
to  haTvi  Nq  by  the  use  of  solution-adaptive  grids,  as  discussed  below.  For  the  implicit  methods,  the 
severe  time-step  stability  limitation  should  disappear.  For  the  explicit  methods,  however,  it  is  inherent 
and  will  remain  a  severe  handicap.  On  the  other  hand,  the  successful  coupling  of  finite-difference  boundary' 
layer  methods,  as  discussed  above,  seems  likely,  and  this  additional  capability  will  undoubtedly  be  very 
attractive  to  future  users. 

Thin-layer  and  full  Reynolds-averaged  Navier-Stokes  methods  -  These  methods  are  considered  together 
because  of  their  similar  numerical  characteristics.  The  main  difference  is  that  the  thin-layer  approxima¬ 
tion  neglects  the  second  derivatives  in  the  streamwise  flow  direction,  thereby  reducing  Ngy  by  about  EOt. 

:t  should  be  noted  that  any  techniques  for  reducing  Wqt  and  Nq  that  are  developed  for  the  implicit  Euler 
methods  above  probably  would  be  applicable  here,  too.  However,  the  greatest  gains  will  come  from  removing 
the  stability  limit  on  at,  which  could  then  be  increased  by  up  to  two  orders  of  magnitude. 

In  practice,  the  principal  reason  for  going  to  the  full  Navier-Stokes  method  would  be  to  capture  some 
separation  phenomenon  for  which  a  finer  grid  in  the  streamwise  direction  would  be  required.  This  is 
reflected  in  the  larger  value  of  Ng,  which  is  the  main  factor  that  makes  the  CPU  time  for  full  Navier- 
Stokes  so  much  greater  than  for  the  thin-layer  approach. 

Figure  19  shows  the  reductions  in  CPU  times  that  could  accrue  from  the  algorithm  improvements  outlined 
above,  plus  the  effect  of  increasing  FLOPS  from  80  million  to  2S0  million.  In  all  cases,  the  number  of 
grid  points  has  been  halved,  and  the  limit  on  has  been  increased  to  0.10. 

It  appears  that  the  nonlinear  potential -flow  methods  will  become  very  economical,  and  that  the  more 
sophisticated  methods  will  no  longer  be  out  of  the  question  *or  specialized  aeroelastic  analyses.  Of 
course.  Fig.  19  does  not  show  the  important  factor  of  how  rapidly  the  improvements  can  be  realized  in  prac¬ 
tice.  For  example,  the  potential  improvements  .qf  the  small -disturbance  approach  are  less,  but  they  can  be 
attained  much  quicker,  than  those  of  the  Euler  and  Navier-Stokes  methods.  Therefore,  we  must  reiterate 
that  the  range  of  ^low  models  described  above  ropresents  widely  Jifferent  levels  of  maturity  and,  hence, 
levels  of  confidence  that  aeroelasticians  are  lively  to  ascribe  to  them  m  practical  applications.  This 
would  seem  to  suggest  that  the  small-disturbance  md  ful  1 -potential  ,.odes  will  remain  much  more  popular  'or 
the  next  ‘ew  years,  despite  the  growth  m  speed  and  memory  if  the  new  supercomputers. 

however,  once  they  have  been  validated,  the  Euler  and  Navier-Stokes  codes  will  help  to  provide  new 
jnderstandi ng  of  and  insights  mto  complev  'low  pnenomena ,  is  well  as  to  generate  lata  bases  which  can  be 


Fig.  19.  Projected  solution  times  for  a  wing  of  moderate  complexity;  75  chords  travel/case. 

Improved  algorithms,  250  mflops. 

used  to  develop  greater  confidence  in  any  of  the  more  approximate  methods.  If  this  new  capability  is 
properly  combined  with  careful  wind  tunnel  experiments,  a  major  surge  in  unsteady  computational  aerodynam¬ 
ics  can  be  expected  in  the  late  1980s,  and  the  projections  listed  in  Table  1  will  turn  out  to  be  pessimistic. 

3.2.2  Grid  generation  improvements 

An  essential  step  in  solving  three-dimensional  aerodynamic  problems  is  the  generation  of  a  suitable 
grid.  This  is  one  of  the  most  rapidly  growing  areas  of  CFO,  and  there  have  been  several  meetings  in  recent 
years  devoted  exclusively  to  grid  generation;  e.g..  Refs.  38  and  39.  The  recent  survey  by  Thompson 

;Ref.  AO)  is  particularly  noteworthy,  and  a  broader  overview  can  be  found  in  Ref.  5. 

Several  contemporary  concepts  for  generating  grids  are  illustrated  in  Fig.  20.  A  significant  advan¬ 
tage  of  some  panel  methods,  and  of  the  transonic  small -disturbance  method,  is  that  the  appropriate  boundary 
conditions  can  be  applied  on  nonaligned  grids;  however,  most  of  the  more  sophisticated  methods  rely  on 
body-conforming  grids.  The  zonal  modeling,  or  multi-block  concept,  appears  especially  attractive  for  com¬ 
plex  configurations. 

The  schemes  for  generating  grids  can  be  classified  into  two  categories,  algebraic  and  differential, 

based  on  the  types  of  equations  that  are  used  to  compute  the  locations  of  the  grid  points.  Although  the 

complexity  may  vary  greatly,  the  algebraic  methods  are  direct  approaches.  The  differential  methods  involve 
the  solution  of  either  the  elliptic,  hyperbolic,  or  parabolic  partial  differential  equations;  this  is  done 
Iteratively  in  the  elliptic  schemes,  or  noniteratively  in  the  hyperbolic  and  parabolic  schemes.  Refer¬ 
ences  5  and  40  include  discussions  on  the  degree  of  control  that  each  method  provides  for  varying  the  mesh 
spacing,  cell  volumes  and  proportions,  and  skewness  of  the  grid  lines.  The  elliptic  and  "nonconformal " 
algebraic  schemes  allow  exceptionally  high-quality  grids  to  be  generated  about  very  general  body  shapes, 
but  they  are  the  most  computationally  intensive  and  expensive.  For  this  reason,  the  more  simple  types  of 
algebraic  schemes  and  hyperbolic  schemes  seem  better  suited  to  the  class  of  unsteady  problems  that  requires 
regenerating  the  grid  at  each  time  step. 

The  concept  of  adapting  the  grid  to  some  feature  of  the  solution,  such  as  clustering  grid  points  in 
regions  of  large  gradients,  has  considerable  potential  for  obtaining  the  maximum  accuracy  with  the  least 
number  of  grid  points.  Therefore,  this  is  an  area  of  active  research  in  the  CFD  cormiuni ty .  As  an  example. 
Fig.  21  shows  a  solution-adapted  grid  and  the  associated  laminar  Ravier-Stokes  results  (Ref.  41)  for  an 
unusual  bluff  body.  The  hyperbolic  scheme  that  was  used  in  adapting  the  inner  grid  to  the  developing  vor- 
ticity  field  required  only  a  small  increase  in  the  CPU  time,  due  to  the  relatively  small  increase  in  Wq-j-; 
but  it  permitted  considerably  better  resolution  of  the  flow  field  to  be  achieved  near  the  corner  of  the 
body.  It  is  work  of  this  type  that  leads  us  to  believe  that  the  grid  point  requirements,  and  hence  the 
CPU  times,  can  be  reduced  by  a  factor  of  two  or  more  without  sacrificing  accuracy,  as  presented  in  Sec¬ 
tion  3.2.1. 

3.2.3  Turbulence  modeling 

The  simulation  of  the  dynamics  of  turbulence  remains  the  foremost  challenge  in  fluid  dynamics  today, 
and  turbulence  modeling  is  probably  the  weakest  link  in  the  chain  of  computational  aerodynamics  technology. 
The  computational  power  available  in  the  foreseeable  future,  in  terms  of  both  speed  and  storage,  will  pre¬ 
clude  adequate  resolution  of  the  broad  range  of  interacting  turbulent  scales  (both  spatially  and  temporally! 
associated  with  most  aerodynamic  flow  fields  at  flight  Reynolds  numbers.  As  a  result,  turbulence  modeling 
has  taken  the  approach  of  single-point  closure  of  the  Reynolds-averaged  Navier-Stokes  eguations,  and  no 
single  turbulence  model  exists  that  can  be  applied  to  i  general  vanetv  of  tlows. 
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Fig.  20.  Grid  generation  concepts  (Ref.  5). 

VORTICITY  CONTOURS  SOLUTION-ADAPTED  INNER  GRID 


As  Shown  in  Table  2  and  Fig.  18,  enormous  computer  resources  are  required  to  solve  time-dependent 
problems  with  finite-difference  simulations  of  the  Reynolds-averaged  Navier-Stokes  equations.  Even  the 
solutions  that  have  been  published  for  steady  flows  have  used  grids  whose  fine  spacing  is  limited  to  the 
single  direction  nearly  normal  to  the  body,  and  hence  fall  within  the  spirit  of  the  thin-layer  approxima¬ 
tion.  This  resulting  computational  process  qualitatively  simulates  separated  flows  and  flows  with  large- 
scale  unsteady  behaviors,  but  the  accuracy  of  such  simulations  is  still  controversial. 

The  recent  survey  papers  on  turbulence  modeling  for  computational  aerodynamics  by  Marvin  (Ref.  42)  and 
Lomax  and  Mehta  (Ref.  43)  and  the  Proceedings  of  the  1980-81  AFOSR-HTTM-Stanford  Conference  on  Complex 
Turbulent  Flows  (Ref.  44)  indicate  the  current  state  of  the  art  in  this  area,  as  well  as  what  will  be 

In 


required  in  the  future.  Since  no  universal  model  exists,  most  researchers  are  now  focusing  their  atten¬ 
tion  on  creating  a  catalog  of  models  based  on  fundamental  Ouilding  block  experiments,  and  most  of  these 
models  are  being  carefully  tested  computationally  to  determine  their  capabilities  and  limitations. 

An  overview  of  turbulence  models,  ranging  from  essentially  no  modeling  at  all  to  the  hypothetical  full 
simulation  of  turbulence,  is  given  in  Table  5.  For  unsteady  flows  of  interest  to  aeroelasticians ,  only  the 
viscous-wedge  representation  of  shock-wave  boundary-layer  interaction,  integral  boundary  layer  (velocity- 
profile  modeling),  and  zero-equation  eddy-viscosity  models  have  been  used  up  to  now.  In  the  larger  domain 
of  steady  flows,  the  primary  research  attention  has  turned  to  multiple-equation  eddy-viscosity  and 


TABLE  5.  TURBULENCE  'iOOELS 


Model 

Physical 
Uenera 1 i  ty 

Numerical 
Compatibi 1 i ty 

Remarks 

Viscous  wedge 

Very  low 

very  high 

Shock-B.L.  interaction 

Integral  B.L. 

Low 

High 

Very  good  when  highly  tuned 

Eddy  viscosi ty 

Zero  eqn 

Low 

High 

Needs  more  tuning 

2-eqn 

Med  1 um 

Low  lu  nigh 

i  20%,  Zt  =  ? 

Reynolds  stress  equations 

Hign 

Low 

3-0  separation? 

Large  eddy 

Very  high 

Low  (?) 

Guidelines  for  above 

Complete  simulation 

Complete 

nth  generation  supercomputers 

Reynolds-stress-equation  models.  Today,  both  experiments  and  specialized  large-eddy  simulation  calcula¬ 
tions  (Ref.  (45)  provide  guidance.  However,  the  current  calculations  of  practical  flows  use  turbulence 
models  that  are  "tuned"  in  conjunction  with  the  numerical  procedure  for  a  specific  class  of  flow  problems. 

As  a  result,  validations  by  means  of  experimental  comparisons  are  mandatory,  and  confidence  in  the  absolute 
values  of  the  numerical  predictions  remains  low. 

In  principle,  the  more  general  models  should  cover  a  wider  range  of  flows  with  less  "tuning,"  but  users 
of  the  large  aerodynamics  codes  may  not  always  feel  that  "bigger  is  better.”  Nevertheless,  it  would  seem 
that  many  of  the  nonlinear  aeroelasticity  problems  will  involve  some  degree  of  flow  separation,  and  in  three 
dimensions.  In  such  cases,  the  two-equation  eddy-viscosity  models  may  turn  out  to  be  the  best  compromise 
between  simplicity  and  generality.  Some  of  the  models  in  this  category  lead  to  stiff  equations,  and  this 
raises  again  the  problem  of  restrictive  values  of  ct.  However,  progress  is  being  made  to  overcome  this 
limitation  (Ref.  46). 

From  this  brief  overview,  it  is  clear  that  turbulence  modeling  will  remain  a  primary  pacing  item  in 
computational  aerodynamics  over  the  next  decade,  for  both  steady  and  unsteady  applications. 

3.2.4  Vortical  flow  model ing 

Whereas  the  treatment  of  shock  waves  in  transonic  flow  was  a  major  focal  point  for  computational  aero¬ 
dynamics  in  the  1970s,  compressible  flow  fields  with  embedded  regions  of  concentrated  vorticity  will  prob¬ 
ably  gain  prominence  in  the  coming  decade.  Figure  22  illustrates  some  typical  cases  where  vortices  interact 
with  components  of  the  vehicle.  Such  nonuni  form! ties  in  the  flow  can  be  expected  to  alter  the  steady  and 
unsteady  loading  and,  hence,  to  have  aeroelastic  consequences. 
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Vortices  sued  from  shd'O  ’eadinq  edges  gave  been  computed  by  a  variety  of  methods;  representative 
examples  can  be  found  m  the  Proceedings  of  the  recent  AGARO  Symposium  on  Aerodynamics  of  Vortical  Type 
flows  in  Three  Dimensions  ;Ref.  A7',  dithin  the  context  of  the  present  emphasis  on  finite-difference  CFO 
methods,  the  direct  approach  ,e.g.,  Rizti  and  Eriksson  (Ref.  A8)  and  Krause  et  al.  (Ref.  49))  may  be 
thought  of  as  vortex  capturing.  That  's,  this  approach  Is  analogous  In  some  sense  to  the  more  common 
sho  .-capturing  methods  that  are  used  m  many  computational  aerodynamics  codes.  As  in  shock  capturing, 
the  details  of  the  actual  phenomenon  are  smeared  over  several  grid  points;  thus,  the  solution  In  that 
region  is  artificially  grid-dependent  and  susceptible  to  the  effects  ot  numerical  dissipation.  In  the 
case  of  concentrated  tip  vortices,  for  example,  the  numerical  dissipation  may  destroy  the  core  structure 
and  Its  large  gradients  ‘aster  than  would  physical  dissipation  (Ref.  50). 

Two  alternate  methods  for  modeling  vortex  flows  are  shown  in  Fig.  23.  By  analogy  with  shock  fitting 
in  transonic  flows,  these  may  be  thought  of  as  vortex  ‘ittinq.  The  upper  part  of  the  figure  portrays  the 
method  of  Caradonna  et  al.  .Re*'  s’  .  whicn  permits  concentrated  vortices  to  be  Introduced  Into  potentlal- 
‘low  formulations.  The  prescr’ped  vonte-  -.jtnod  shown  in  the  lower  half  of  the  figure  is  due  to  Steinhoff 
,Ref.  il':,  and  it  nas  also  been  jsed  iuvcess*u 1 1  v  by  Snn'wasan  and  McCroskey  (Ref.  53)  for  potential, 
Euler,  and  thin-layer  Navier-Stoxes  analyses  ;*  ..nsteadv  1 1 rfoi  1 -vortex  interactions.  In  these  applica¬ 
tions,  the  structure  of  the  vortex  's  prescnped,  but  ^ts  oath  m  space  develops  as  part  of  the  solution. 
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L'-eai  pane’  me'nogs  e  j.  .  ''asxew  Ref.  54  ;  and  vor*ex-fi  lament  methods  (e.g.,  Rehbach  (Ref.  55)  and 
^eonarg  Re*',  to  nave  been  jsed  t.i  treat  vortex- dominated  ‘lows.  However,  these  techniques  have  not  been 
applied  tj  compressible  ‘' ;ws  jP  to  now.  'he  cacnOination  of  shock  waves  and  regions  of  concentrated  vor- 
ticity  represents  a  challenging  put  ‘ruit‘u!  area  of  research  over  the  next  few  years. 

l.^.b  Aerodynami c. structura 1  coupling 

It  's  sometimes  li‘ficult  ‘or  tamoutationa '  aerodynamici sts  to  realize  that  their  Impressive  results 
are  merely  the  ‘orcinq  functions  ‘.or  complex  lynamic  systems,  and  that  aeroelastic  analyses  are  often  domi¬ 
nated  by  structural  cons i dera t i ons  and  calculat'on  methods  'hat  are  guite  alien  to  the  world  of  CFD.  And 
as  we  nave  noted  previously,  few  specialists  are  highly  skilled  both  ’n  computational  aerodynamics  and 
structural  dynamics. 

'he  traditional  method  )f  predicting  ‘’utter  boundaries  has  been  to  determine,  first,  the  unsteady 
iirioads  via  aerodynamic  influence  :oef ‘’cents .  'hese  could  be  computed  from  linear  theory  as  functions 
;f  the  geometry  of  the  wing,  the  fach  number,  the  reduced  ‘requency.  and  the  simple,  uncoupled  motion  of 
the  wing,  'his  relatively  simo’e  approach  ;an  still  be  used  -n  the  transonic  regime,  if  the  aerodynamics 
can  be  considered  as  time  linearized,  and  relatively  modest  comoutinq  power  is  required. 

If  the  unsteady  components  of  *he  airloads  depend  nonlinearly  on  the  structural  motion,  then  some  sort 
of  direct  coupling  is  required  and  time-accurate  calculations  must  be  performed,  as  indicated  In  Fig.  24. 

In  this  case,  the  computer-resource  requirements  r’$e  sharply  Two-dimensional  calculations  of  this  type, 
in  which  the  itrijctural  equations  of  motion  and  the  'ransonic  sma  1 1 -di sturbance  equations  were  integrated 
simultaneously,  have  been  reported  Oy  Huruswamy  and  ‘anq  'Ref.  57',  ‘or  example. 

If  the  series  or  sequential  couplinq  ’echnique  it  ‘’g.  24  's  used,  several  iteration  cycles  for  each 
‘''■ght  condition  may  be  required  to  determine  the  ‘'utter  boundaries.  The  results  can  be  obtained  faster 
if  vector  processinq  or  pipelining  techniques  are  used,  but  the  solution  methodologies  for  the  structural 
and  aerodynamic  parts  can  be  dif‘rrent.  On  the  other  hand.  i‘  parallel  coupling  were  used  within  each 
time  step,  the  number  if  iterations  would  probably  be  'ess.  However,  the  two  parts  of  the  program  would 
have  to  be  totallv  lomoatible  and  wou’i  have  to  be  coded  very  arefully  if  the  computer  resources  were  to 
be  jtii’zed  effect 've’y. 

A*  'east  ‘or  'he  next  *ew  vimx.,  .  *  ’ ,' i  te- 1 1  “erence  methods  ire  'iKelv  to  remain  the  primary  tool  fpr 
analyzing  ,,omple<  jnsteadv  ‘’ow  *’e’.js  ing  * .  r  ■  tp-ei  emenr  *r,f  .ina'yzina  complex  structures,  and 
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24.  Coupl’nq  of  derodynamic  and  structural -dynamic  analyses. 

the  compatTbi 1 i ty  of  these  two  iooroaches  Tiay  be  a  oroblem.  A  possible  strategy  that  would  allow  both 
oomputational  aerodynamicists  ana  structural  dynamicists  to  concentrate  fnore  in  their  own  specialties  is 
shown  'n.  Fig.  25.  "his  aoproacn  ^ould  use  multiprocessor  supercomputers;  one  processor  solves  the  aero¬ 
dynamic  equations  by  whatever  ^thod  ’S  the  most  e*^icient  and  appropriate,  while  the  other  processor  is 
working  :n  the  structural  equations,  and  they  would  (perhaps)  share  the  large,  comnon  memory  of  the  central 
computer  facility.  >  any  case,  the  jeneral  ’SSue  of  efficiently  coupling  the  aerodynamic  and  structural 
parts  of  the  *’^utter  oroplem  wil’  r-emam  a  ma:or  challenge  for  the  foreseeable  future. 
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F^g.  25.  Ae'’Odynam7 '••jctur  j  1  :/uo’inq  jsinq  multiprocessor  supercomputers  . 


3.2.5  Display  and  ore-  inq  post-urocessin;}  3ata 

Simply  defining  the  leome'.r/  >*  j  ■.omple''  rr-.raft  int’ jurat'on  can  be  a  formidable  task,  as  attested 
to  by  a  '■ecent  symposium  jn  this  .uD.ect  ^ef  Sd : .  addition,  the  solution  of  three-dimensional 


aeroelastic  problems  Ti/ol.es  tne  nanaqement  of  encrmous  amounts  of  data,  especially  output  data.  This  can 
significantly  affect  the  productivity  of  the  analyst  who  nas  to  consume  that  data,  make  decisions  about  it, 
and  formulate  new  avenues  pf  approach.  Therefore,  the  need  to  digest  efficiently  this  quantity  of  data 
makes  the  development  of  optimal  pre-  and  .oost-data  processing  procedures  absolutely  essential.  Pre-, 
intermediate-,  and  oost-processmg  of  bulk  data  can  only  oe  done  effectively  using  high-resolution,  high- 
throughput  computer  graphics  devices,  ’'hus,  the  efficient  use  of  on-site  supercomputers  will  necessitate 
networking  them  with  peripheral  minicomputers  whicn  are  linked  with  sophisticated  interactive  graphics  work 
stations,  as  indicated  in  Pig.  2b. 


<  c 

-t:i.  .6.  Tatj  cocessino  and  analysis. 

'leneration  pf  three-dimer  lonal.  time-vary: n-;  ’at.i  coses  /et  another  significant  problem  for  the  aero- 
elast'cian,  namely,  now  to  .Piso'.iy  tne  data  ■'leanmcgfu' iv  ‘‘or  potimum  understanding  and  analysis.  Modern 
graphics  displays  are  improving  ,■  a  rapid  pace,  and  two-  and  three-dimensional,  computer-generated  color 
pnotograohs  can  be  used  to  ireat  idvantage  Pefs  t9  and  62).  Figure  27  is  reoresentative  of  current  capa¬ 
bilities.  The  wind  tunnel  •'Ode'  ■•‘t  .'a  ■  's  ^'rs:  '-..r^resented  by  tne  wire  frame  with  hidden  lines 

removed  F'i.  27b''.  .  th"  haoeo  .i.ir"!  ,  •  !■  eteg  sur*'ace  t  eoresentation  [Fig.  27c)  is  much  more 

vaiuaole  m  /isua’-jmg  the  node’  ang  i-  tpottirq  errors  'n  the  lurface  representation.  Finally,  the 
,o1or-shadi  ng  rporesenta  1 1  sn  pf  the  ■'.tmpu'ed  pressure  ii  ■■.  tr-pution  on  the  surface  'Fig,  27d)  is  an  impor¬ 
tant  aia  'n  evamininip  the  , elution. 

jraphicai  display  iev'ces  mg  supporting  .oftware  ire  now  available  to  generate  three-dimensional 
color  movies,  'n  both  ing  /ideo.  These  car  be  invaluable  'or  visualizing  an  evolving  sequence  of 

aerodynamic  events,  ,ucn  as  ''.jtter.  and  ‘.jr  effect'vely  presenting  and  describing  these  events  to  the 
tecnnica’  .c.pnmunity. 

, n  a ‘d d  1 1 1 'J n  ,  ma i.^'r^s  are  al>Ci  j/diidult?  ‘  jr  /  Droduc  ’  ng  hardcopy  directly  from  the 

graphics  device  ’•or  i  copy  instant  viewing  j  color  -ncvie  from  video  disk  or  tape. 

However,  the  prcduct’-.n  o*  uch  'cior  -lov’es  ,  pen'i  t  ve .  Although  *'ui  1  color,  ’‘.hree-dimensional  hidden 
ime  TKiv’es  wil'  ’■eivi'-e  , '  ant  •  .rTpu'er  -♦•ne.  t  'idy  oe  ;  ns  igni  ^'uant  when  compared  to  the  computer 
time  reoui’^ed  tc  ;enf'*dte  the  lata  t  >#n;cn  novie  j  lade.  furthermore ,  the  oenefits  to  be  gained 

,uch  \  nqv’r^  ■  "b  .rodut'  on  .c.t, 

j.J.’  Arti^it.'a’  n  te  '  '  ■  f' .  •<  now’ ►'doe -oased  ./S'o'tis 

>.  the  pr»ced’ng  ue  'Mn-,.  '^dvn  oi  .a.jsued  *hp  .-lifter  Hardware  and  software>  the  complex  algorithms* 
the  gr'js.  the  .'dh/s'ci’  ^ode’^n.;,  jerodvramv  tructur,:  .-uol’ng,  and  the  data  processing  that  will  be 
•nvolved  ’n  aeroelastic  sinulat^cns  *0  ’he  -hese  iimulations  will  require  vast  expertise 

and  enormous  resoijrres  in  ".erms  t<otn  ^umar  ^or^  -'"-puter  rapacity.  The  aeroelastician  will  need 

highly  developed  skills  'c  *:he  d-uC’pl’ne,  -*  jn«',s:s,  lerodynamics*  structural  dynamics,  com¬ 
pute’*  jCience,  and  time  Tanageme’^t.  '/  ‘hf:.  rv  •  .ncepts  of  irtificial  intelligence  (AI)  can 

be  jopned,  especially  the  concept  it  •M-.-Dr;:  .;y .  .ugqested  in  Refs  61  and  62. 

expert  systems  are  knowledge-based  SvSrpu'’  ’ndt  perform  specialized  tasks  at,  or  perhaps 

beyond,  *he  level  of  a  human  expert,  'hi  v  ”1:0  '»'ve'  *  •'nant  f?  is  a  ’'esult  of  "domain-specific”  knowl¬ 
edge  and  strategies,  expert  j/Stem$  are  -'i.'  .vstems  are  d  i  s  ti  ngui  shed  from  other 

artificial  intelligence  .irnqrams,  and  rmpu-er  .  rograr^ ,  'f  jf'-t-'-j',  by  their  ‘bilitv  to  reason  about  their 
own  processes  of  inference,  md  "o  “'ur'M  .n  p.;;!  v’lt  •  f '-*;.irrr/;  'n. ■,',»>  iirncesses.  These  distinguishing 
v:  jractenstics  are  made  possibi*-  ov  'ht*  jnJe* :  -rr  f' ;  *L;rt‘  .amnn  •'o  most  expert  systems.  As  indi- 
ated  m  Fig.  28,  the  two  major  ■  mptinent  >  arr-  ,  ,.wl 1  :»•  ‘MbH.  ^ni.  h  consists  of  domain-dependent  facts, 
nijles.  jnd  heuristics,  and  a  .eDdrjte  •  rro’ cr- .  .vctem  ’0  proceed  efficiently 


blatters  tnat  iet  eN2t>r-ts  icart  ’■rom  aegmners,  are  symoolii;,  infe''ential ,  and  rooted  in  experiential 
knowledge'  Ret.  63).  As  indi>.atea  in  rvg.  38,  three-dimensional  grid  generation,  nonlinear  flutter 
analysis,  and  data  processing  are  three  areas  that  seem  guite  compatible  with  these  characteristics. 

The  knowledge  base  for  such  aoplications  would  consist  of  facts  such  as  flow  solvers,  structures 
solvers,  turculence  models,  vortex  mode’s,  and  grid-generation  schemes)  and  heuristics  ii.e.,  experience 
and  good  judgment  concerning  flutter  characteristics  of  aircraft).  The  inference  procedures,  wh  .  cannot 
be  defined  precisely  at  tnis  point,  would  process  the  information  that  is  calculated  by  the  various  stages 
of  the  program  and  would  make  decisions  and  draw  conclusions  based  on  the  successive  accumulation  of  new 
facts,  starting  with  previous  solutions.  An  important  feature  of  the  expert  system  is  its  ability  to 
process  information  from  many  'outside'  experts,  as  indicated  m  fig.  38. 

It  .must  be  empnasized  that  constructing  an  expert  system  requires  a  suostantial  investment  of  time 
and  manpower.  Estimates  va  ■>  anywnere  n  ’ess  than  a  year  -'or  simple  systems  in  a  friendly  environment 

with  existing  tools,  to  15  years  "'or  comp  e-  systems  m  lamanoinq  environments  where  fundamental  research 
and  development  is  '■eouired.  furthermore,  the  'evel  of  performance  of  various  systems  vanes  greatly. 

Some  of  the  systems  tnat  already  nave  been  puilt  now  routinely  solve  practical  oroplems,  while  others  have 
never  made  't  beyond  the  research  stage  A  necessary,  but  not  sufficient,  condition  for  success  seems  to 
be  that  experienced  researchers  ana  mature  tecnnologies  must  be  used  ♦or  builoing  almost  any  expert  system. 
In  any  case,  the  potential  payoff  'rnm  the  ,se  or  expert  systems  ih  the  numer'cal  aeroelastic  simulation 
process  ’S  worthy  of  jttention,  ana  't  warrants  the  allocation  jf  resources  as  an  investment  in  the  future. 


d.  Eij'MMARY  AND  EDNCUSIONS 

Eomou tat  1 una 1  aerodynanics  's  now  a  well-established  tool  'hat  is  used  throughout  the  world  in  the 
design  ana  analysis  jf  modern  flight  vehicles,  'he  past  decade  has  -.een  impress. ve  growth  in  the  capabili¬ 
ties  ''or  oredictinq  itatic  performance  and  airloads  aver  a  wide  '.-.qe  of  conditions,  as  evidenced  in  the 
examoles  given  m  Section  3.  "urtnermore,  this  growth  seems  to  be  sustaining  itself,  both  in  terms  of  com¬ 
puter  speed  and  memory  and  in  the  numerical  methods.  Exceot  for  tne  I’mitations  of  turbulence  modeling, 
the  showcase  problems  .gf  ’97A  can  be  solved  routinely  today,  br-oqress  in  the  more  specialized  areas  of 
unsteady  aerodynamics  ana  flutter  analysis  has  tended  to  lag  the  jevelooments  m  steaoy  aerodynamics;  but 
the  .dace  is  quickening,  and  the  potential  'or  further  improvements  now  seems  to  be  even  greater, 

'he  primary  factors  that  limit  computational  aerodynamics  'n  its  aoDlications  to  aeroelastic  analyses 
are  tne  costs  of  performing  a  large  number  of  nonlinear  calculations,  "id  the  validity  of  the  numerical 
simulations  The  validity  is  essentially  determined  bv  the  turbulence  modeling,  by  the  ability  of  the  grid 
to  resolve  the  relevant  details  of  the  'low  field,  and  by  the  accuracy  of  the  finite-difference  solutions 
■h  reoresenting  the  physical  flow.  The  cost  of  the  calculations  is  determined  by  the  computer  hardware  and 
software;  by  the  manpower  reouired  to  implement  the  codes  and  to  digest  the  "esults.  and,  in  typical 
unsteady  ca '  cu  i  a  tions ,  by  the  stability  restrictions  on  the  time  steo  that  ,;an  be  used  for  the  ■more  sophis- 
ticateJ  methods  df  analysis. 

'ortunatel y .  tne  trends  are  highly  favorable  ♦or  most  of  these  limning  ‘actors.  The  area  of  turbu¬ 
lence  mode’ing  is  proDaply  tne  one  with  the  least  optimism,  although  the  manpower  limitations  are  of  con¬ 
cern.  'ho  ana’ysis  in  Section  3  indicates  that  a  complete,  and  perhaps  guant i ta ti ve .  simu'ation  of  the 
transonic  viscous  ‘low  oven  complex  configurations  will  become  possible  within  the  next  decade.  This  goal 
wii'  be  ach’eved,  however,  only  with  the  aid  of  hi gh-dual i ty  onysica’  experiments.  That  is,  detailed 
exper-rents  wi’l  nave  to  p’ay  crucial  'dies  in  improving  the  turbulence  and  vor'ex  modeling  and  in  guiding 
and  validating  the  numerical  simulations,  whatever  their  levels  of  complexity. 

However,  the  cost  to  calculate  hundreds  of  combinations  of  ‘low  parameters,  structural  frequencies  and 
mode  .ihaoes,  and  w-g-store  configurations  will  be  high  enough  so  that  aeroelasticians  will  continue  to 
want  ’’ono  aporoximate  'engineering"  methods.  The  develooment  of  better,  and  less  costly,  more  approximate 
techniques  wil'  be  greatly  enhanced  by  an  intelligent  combination  of  the  large-scale  numerical  simulations 
and  wind-tunnel  exoeriments ,  Finally,  the  expert-system  concept  of  artificial  intelligence  could  possibly 
hasten  tne  achievement  of  the  aeroelasticians'  goals. 
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